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Abstract 

The behavior of the diameter of the coexistence curve in terms of the en- 
tropy and the corresponding diameter are investigated. It is shown that the 
diameter of the coexistence curve in term of the entropy is sensitive to the 
change in the character of the rotational motion of the molecule in liquid 
phase which is governed by the short range correlations. The model of the 
compressible effective volume is proposed to describe the phase coexistence 
both in terms of the density and the entropy. 

1. Introduction 

Peculiarities of the asymmetry for the vapor-liquid binodal are the im- 
portant source of information about different factors influencing its form. 
It is well known that the asymmetry of binodal is absent for the simplest 
system known as the lattice gas. From physical point of view it is caused 
by 1) the absence of the thermal expansion effects (the lattice distance is 
constant) and 2) the invariance of the energy of system about the trans- 
formation: particles ^ holes. The latter circumstance is evidently violated 
by 1) the appearance of the soft repulsive core effects and the long-distance 
attractive interaction between molecules, which are caused by the dispersive 
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and electric multipole-multipole interactions, as well as 2) the anisotropic 
character of interactions for non-spherical molecules. These effects lead to 
the essential difference in the coefficients of the thermal expansion for liquid 
and vapor coexisting phases and violate the particle-hole symmetry. One of 
the simplest characteristics of the asymmetry of the vapor-liquid binodal is 
its diameter. 

If the vapor-liquid binodal is described in the terms temperature-density, 
its diameter is determined as: 

where rij, i = I, v, are the densities of the hquid and the vapor phases cor- 
respondingly and He is the critical density. The number density (or specific 
volume) is the natural order parameter for the liquid-vapor critical point 
(CP). In its vicinity the asymmetry of the density diameter is mainly deter- 
mined by the long-range fluctuations of the order parameter [|T|. Far away 
from the CP the asymmetry of the diameter of binodal is mainly caused by 
the hard core effects. This statement follows immediately from the van der 
Waals (vdW) equation of state, which gives us the rehable qualitative infor- 
mation about the role of main factors in the formation of the EoS. However, 
the description of the binodal in terms of the density is rather incomplete. 
The density is expressed via the unary correlation function. Obviously, at 
such level the rotational degrees of freedom as well as the short range correla- 
tions due to nonspherical shapes of the hard core do not manifest themselves 
in a clear-cut way. 

The information about the role these factors can be obtained from the 
analysis of the entropy diameter: 

S, = ^^-S. (2) 

In general, it includes the part changing similarly to density and the irre- 
ducible part, which is caused by the correlation effects of all orders. Due to 
this the analysis of the asymmetry of the binodal in terms of the temperature- 
entropy is more informative in comparison with its temperature-density ver- 
sion. 

This circumstance becomes evident from the comparison of the temper- 
ature dependencies of nd{t) and Sdit) for several classes of substances: 1) 
the noble fluids; 2) the simple molecular fluids with the inversion center for 
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their molecules; 3) the simple molecular fluids without the inversion center 
for their molecules and 4) the high- molecular fluids. The behavior of nd{t) 
and Sdit) for these of substances is presented in the Fig. |1| and Fig. |^. From 
the Fig. it follows that the density diameter is the monotonous function 
of the temperature at least outside the fluctuational region. The tempera- 
ture dependence for all enumerated substances is close to linear in the wide 
temperature interval, in connection with this it is called as the law of the 
rectilinear diameter p[. In contrast to it the behavior of the entropy di- 
ameter is essentially non-monotonous. Such a discussion testifies about the 
important role of the caloric part of the entropy, which is not reduced to the 
density effect. On the other hand it gives the additional information about 
the molecular motion, in the first place, the rotational motion of molecules, 
which directly contributes to the heat capacity and the entropy. We see that 
the near critical behavior of the entropy diameter is also stronger expressed. 
In IQ the analysis of the behavior of Sd was used to estimate the degrees of 
the dimerization in the noble fluids near their critical points. 

In this paper we investigate the nature of the asymmetry of the binodal 
in the temperature-entropy terms for molecular substances. The entropy di- 
ameter Sd is used as the measure of the asymmetry. We will show that the 
surprising non-monotonous behavior of the entropy diameter is mainly ex- 
plained by the rotational motion of molecules. The interconnection between 
the rotation of molecules and the excluded molecular volume is established. 
The model equation of state (EoS), taking into account the dependence of the 
effective excluded volume on pressure, is proposed. Due to this the behavior 
of Sd as well as the form of the binodal in terms of the temperature-density 
in the wide temperature interval beyond the fluctuation region near the CP 
are reproduced self-consistently. 

2. The qualitative analysis of the non-monotonous behavior of the 
entropy diameter 

In this section we discuss the main physical mechanisms generating the 
non-monotonous temperature dependence for the diameter of the entropy. 
For this purpose we start from the general expression for the specific (per 
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particle) entropy of arbitrary system: 

5 = 5, + cJn^ + /(n)-/K). (3) 

J- c 

Here Sc is the value of the specific entropy at the CP, c„ is the dimensionless 
specific heat {ks = 1), f{n) is the function which depends on the density 
only. 

From (^) it follows that the diameter of the binodal in terms of 
temperature-entropy is determined by the expression: 

S^ = l {c^:Ht) + c^:\t)) In I + 1 (/(n,) + /K)) - /K) , (4) 

where Cv\ i = l,v, are the specific heat capacities on the liquid and vapor 
branches of the binodal. 

Since Cv\t),c^j\t) > 0, the contribution of the first (caloric) term to the 
diameter of entropy enters with negative sign for all temperatures T < T^. 

The numerical values of the combination | (^Cv\t) + c^\t)^ depend on the 

number of the degrees of freedom excited in the thermal motion. For the 
noble gases, having only translational degrees of freedom, 

3<l{cl^\t) + c^:\t))<Q. (5) 

The lower limit corresponds to the case of non-interacting particles, the 
upper - to the strong interaction between them. For the substances with 
non-spherical particles the contributions of the rotational degrees of freedom 
should be taken into account. In this case, 

Q<l{c';!\t) + c^:\t))< 12. (6) 

These inequalities are illustrated by Fig. ^. The vibration degrees of free- 
dom for many substances at T < Tc and, in particular, for all substances 
enumerated in Fig. |l] and Fig. ^ can be ignored. 

The density contribution f[n) to the entropy obviously increases if the 
density decreases. Far away from the critical point f{ni) < (<<)/(n„). The 
resulting sign of Sd depends on the relation between the caloric and density 
contributions. To get the qualitative representation about their competition 
let us consider the ideal gas approximation for /(n). 
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2.1. Ideal gas contribution 

In the ideal gas approximation we should put 

/(n) = -ln-. (7) 

This approximation correctly reproduces the tendency in the change of the 
entropy, caused by the density contribution, for the vapor and liquid branches 
of the binodal. Therefore, we hope that the expression 



sf'W = 1 



(c<"(«) + cW(())lnI-lu"'"" 



Tc n 



2 



(8) 



will correctly reproduce the qualitative peculiarities in the behavior of Sd{t). 
Here we put additionally 

and 

cl'^t) = ci'^Hi) + Ac„ (9) 

where i = itr + ir is the sum of the translational and rotational degrees of 
freedom. Of course the value i is determined only in low density vapor phase 
where the interaction contribution to the specific heat is negligible. With 
the increasing density this contribution increases (see Fig. |^). In accordance 
with Fig. 1^, we put 

{1, for Ar, 
2, for H2S, (10) 
5, for H2O. 

Substituting the estimates (y) and ([TD|) for the specific heat to (^ and the 
experimental values of ni and we find that 1) at z = 3, that is characteristic 
for argon-like liquids, the diameter of entropy is a monotonous function of 
temperature and 2) the non-monotonous behavior of Sd{t) is observed for 
i > 5. These results are in full qualitative agreement with experimental 
results presented in Fig. 0. 

Here it is relevant to pay attention on the position of the roots of the 
equation S^\t) = 0. It has two different roots: ti and t^, ti < t^. One of 
them, tfi, coincides with the CP: = tc = 1. Another root ti is close to 
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Figure 3: The specific heat along the gas branch of the binodal for the number of molecular 
liquids: 1 - Ar, 2 - 02,3 - H2S,4 - NH3, 5 - HjO, 6 - SO2. The dashed lines correspond 
to i = 3, i = 5 and i — 6 respectively. 

the corresponding experimental values. For example, for O2 Eqs. and 
with 2 = 5/2 and Ac^ = 3 gives: 



At the same time, the equation S^^"^^ = has three roots: ti,tu and t/i = 1. 
The roots th for equations Sj^^ = and S^i"^^^ = are trivial, since the 
diameter of the entropy vanishes at the CP because of its definition. For O2 
the intermediate root tu = 0.88. The appearance of this root for the equation 
Sj^^"^^ = is connected with the influence of the critical fluctuations. 

2.2. The manifestation of critical fluctuations 

The critical fluctuations manifest themselves in two relations: 1) they 



shift the position of the CP {Tr^> ^ = Tr^' - ATfi, ATfi > 0), 



determined in the framework of some mean-field approximation and 2) they 
renormalize the value of the entropy at the CP: Sc"^^^ Sc = S^^^ — 
ASfi, ASfi > 0. In consequence of this at the comparison of S^"^"^^ and 
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Sf\ which can be considered as the simplest mean-field approximation, we 
should take into account that the normalized temperatures are different. We 
have 

sf" ( I) « sr i ^.„J \ « sr' { -LI + AS,,. (12) 



TJ ' \t1"'"-AT,iI " \Tt'> 



where 



is the fiuctuation contribution to the entropy at the CP. 

Let us consider the application of Eqs. (0),(0) to the description of the 
diameter of the entropy outside the fiuctuational region for noble gases. For 
the comparison of (t) corresponding to the experimental data with S^j^^ 
it is necessary to shift the latter curve 1) upward on the distance ASfi and 
2) along the temperature axis on —AT/-/. These operations will be carried 
out everywhere below. As a result, the height of the plateaus (see Fig. ^a) is 
determined by the second term in Eq. (p^. At the same time it means that 
the caloric contribution into 5^™''^^ is practically compensated by the density 
dependent term in the temperature interval 0.65 < T/T^ < 0.95. 

2.3. The diameter of the entropy for the van der Waals model of the liquid 
and vapor branches of the hinodal 

In this subsection we consider the temperature dependence of the diam- 
eter of the entropy if the binodal is modeled by the van der Waals equation 
of state. 

P=^^-an\ (14) 



The values of ni and n^, which follows from ( |I^ ) should be substituted in 
(P). Near the critical point we have 



l±5o|r|^/' + v4|r| + ..., (15) 



where r = From ([T^) it is easy to get: 

Bo = 2, ^ = 0.4. (16) 
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Substituting Eq. (|5]) into Eq. with (0) in the whole temperature interval 
for hquid states we obtain: 

Sf\T) = (cf ) + Iac.) In I - In (1 - A|r|) , (17) 

where A = Bq—2 A. For the van der Waals EoS A = 3.2. The same procedure 
apphed to the EoS with the Carnahan-Starhng contribution: 

p = nT an (18) 

{i-vr 

where t] = nb/A, leads to (|T7|) with A ~ 3.7. The comparison of the exper- 
imental data for H2S and the results obtained from Eq. ([T7| ) is on Fig. ^. 
The fitted value of the parameter A ~ 1.8 is almost as twice as smaller than 
the corresponding values which have been obtained above within simplified 
model d^) and (|15]). This simple treatment allows to conclude that the non- 
monotonic dependence of Sd can be get qualitatively as the consequence of 
the interplay between two opposite tendencies: increase of the configurational 
chaos of the gaseous phase and the entropy decrease of the liquid phase with 
lowering the temperature. In Section ^ we propose the model EoS describes 
the behavior of the entropy diameter more adequately (see Fig. |^ below). 

2.4- The excess entropy 

In connection with said above it is appropriate to use the excess entropy: 

5(-) = 5(T,n)-^ln^ + ln-, (19) 

2 Tr Tin 
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in which the ideal gas contribution is subtracted. In accordance with its 
definition S^^^^ includes the fiuctuation contribution in the vicinity of the 
critical point and the effects caused by the interactions between translational 
and rotational degrees of freedom. It means that S^^^^ is expected to be 
greater in absolute value for the liquid phase in comparison with the vapor 
phase. This statement is naturally confirmed by Fig. |. In correspondence 
with this we redefine also the diameter of the entropy 

5, -> = 1 5;^") + Si'^^ ) - Si'^^ . (20) 

It is behavior is presented in Fig. |^. The latter is monotonic excluding 
small interval near the critical point where the fiuctuation effects lead to 
the appearance of singularities 0. The slope of S^^^ is the source of the 
important information about the combination of coefficients describing the 
attraction and repulsion between molecules. 

Here we want to note the following peculiarities of the temperature de- 
pendence of S'f^^ for different liquids. We see that the characteristic values 
of slopes of S^^^"^ for noble gases are essentially smaller in comparison with 
those for molecular liquids. Besides, the values of S^^^ for different noble 
gases are very close that is in good agreement with the Principle of Corre- 
sponding States. The difference between noble gases and molecular liquids 
has also natural explanation: the intermolecular interactions in molecular 
liquids are stronger than in noble gases. 
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Figure 6: The behavior of the diameter of the excess entropy. 
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It is also essential that the slopes of Sf^ for molecular liquids of different 
types are close to each other i.e. the molecular sizes, electrostatic multipole 
moments and H-bonds are responsible only for fine details. Such a situation 
reflects the fact that the behavior of S^^^ is determined by the averaged 
intermolecular potential 0, After averaging on the rotational degrees of 
freedom the effective intermolecular potential reduces to argon-like one which 
in general is the state dependent quantity ^, 0. The noticeable deviations 
from the predictions made with the help of averaged potential are observed 
only in the narrow vicinity of the critical point and in supercooled states 
of liquids. As an example, for water the contribution of the orientational 
correlations can reach several percents [11 . 

3. Effective volume model 

In this Section we propose the model of the effective volume which allows 
to include on the phenomenological level the effects of the nonsphericity of 
the particles. As shown in ||T2Hl4[ the form of EoS ( ^2]) with a and b as 
constants is adequate in the vicinity of the triple point. But they fail to 
reproduce the behavior for Sd even for noble gases like Ar, Kr and Xe where 
the interaction is obviously spherically symmetrical. For the molecular fluids 
with anizotropic hard core potentials like and using vdW-like EoS 
which is based on the effective spherically symmetrical interaction can be 
justified |^. Therefore we can use the simplest geometrical representation of 
the shape of the excluded effective volume as a sphere. 

Nevertheless such a fine property as the behavior of Sd was not studied. 
In additions the specific heat for the model EoS's of type (^) with b = const 
does not depend on the density. As a result the diameter of the entropy 
as well as the behavior of the specific heat along the coexistence curve is 
not reproduced even qualitatively. To overcome this deficiency one should 
consider the parameter b as the function of the thermodynamic state. 

In this Section we propose simple model which account for the rotational 
degrees of freedom of the non spherical molecules through making the pa- 
rameter b state dependent quantity. The correspondingly modified EoS of 
Carnahan-Starling is proposed. The first consequence of such a modifica- 
tion is more correct reproducibility of the liquid branch of the binodal. The 
second is the correct description of the entropy diameter and satisfactory 
agreement with the specific heat data outside the fluctuational region. Note 
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that van der Waals himself and Kamerhng-Onnes did attempts to make the 
vdW parameter b pressure dependent [p!^-|T7 . 



From what has been said above it is clear that one can model the influence 
of the density on the rotation as the restriction of the angular conflguration 
space available for the molecule in a cell formed by its neighbors. The change 
in the available space for free rotation of the molecule can be described by the 
dependence of the volume parameter b on the density and the temperature. 
In such an approach the change of the excess entropy of the liquid phase with 
the density (see Fig. ^ can be reproduced. 

Basing on the physical interpretation of b as the available free volume we 
propose the following model: 



1+70; 



X 



piid) l^(vdW) ^ pM^^rp^ (21) 



where p^c'^^^ is the vdW critical pressure and 7 is the adjustable parame- 
ter. From the physical point of view the model (pl|) states that the volume 
available for a particle decreases if the pressure increases: 

f = fO (1 - IPid/PQ + •••)• 

This qualitatively conforms with the pressure dependence of the free volume 
which can be obtained from the molar refractivity data JTSf . The conception 



of the effective free volume in the liquid state theory is widely used. A lot of 
modifications for the repulsive pressure contribution which use temperature 
dependent effective core diameter a{T) have been proposed |IS| -PT[]. Usually, 
the P — V — T data are used for the determination and the interpretation of 
the effective free volume and in particular the soft core diameter cr(T). The 
analysis of the influence of the core softness on the entropy and the speciflc 



heat are more scarce 



To justify the model ( pID we use the generalized vdW EoS approach 



23| , |2^ . This approach allows to deflne the free volume in rigorous manner 
for the molecular systems with the spherically invariant interactions. It is 
based on the fact that for a wide class of the potentials the EoS for simple 



molecular fluids is PH]: 



p = p+(n,T)+p_(n,T), (22) 

where p+ is the pressure contribution due to hard core repulsive interactions 
and p_ is the contribution of the attractive long range part of the potential. 
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Two contributions in (^) can be represented in the vdW-like form 

nT , . 

1 — on 

p_ = -av? . (24) 

The vdW coefficients a(n, T), 6(n, T) are in general the functions of the ther- 
modynamic state PB|. The function b{n,T) is determined by the contact 
value of the cavity function: 



1 + bony{a,n,T) 

in accordance with the basic representation for p5| : 

p+ = Pidil + bony{a,n,T)) , (26) 

where Bq = = 4f o and a is the diameter of the hard core, y{r, n, T) = 
e^'"^'^^ g{r,n,T) - the cavity function, g(r,n,T) - pair distribution function. 
Here Pid = nT is the ideal gas pressure. E.g., for the Carnahan-Starling (CS) 
model of p+ the function b is determined as: 

t = ifi ^ 



1 - , ,1 , ,3 1 ^ V = nbo/i. (27) 



In [^] it is shown that in the high density regime n > n*{T), when the 
parameter a{n,T) becomes negative {a{n*{T),T) = 0), the approximate EoS 
has the form 

pil-bfu) = nT (28) 
and the free volume of a molecule bj can be written as: 

(29) 

p J 

Comparison ( PT| ) with ( PP] ) justifies the model of compressible volume as the 
simple approximation for (pQ]). 



Note that the structure of (|26|) is determined by the form of the EoS (p2| 
due to separation of the hard core effect. Indeed, the ideal gas term in (^ 



is the zeroth approximation if 6 — )■ 0. But as noted in |2^ the parameter b 
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(a) (b) 

Figure 7: The comparison of the resuhs for the binodal and the diameter of the entropy 
for H2S calculated within proposed model of compressible volume Eq. with 7 w 0.03 
(solid lines) and the common assumption h = const ,7 = (dashed lines). The points are 
the experimental data. 



is rather the excluded volume than the size of the molecular core although 
there is correlation between these quantities [|18|. Therefore we can extend 
([26| ) by introducing the compressibility of the excluded volume h into any 
model for in particular for CS approximation (^). The corresponding 
pressure term has the form: 



ACS) 



n dh 



(30) 



where p^^^^ is the repulsive pressure term in CS approximation, which now 
corresponds to the limit of incompressible effective volume — ?■ in 
analogy with (|^). In accordance with Eq. (|30|) and the identity p = — |y 
we arrive to the following free energy consistent with Eq. (pO]): 



^bn( 4 — bn) , . 

F = F,, + T ^ ^ -an, (31) 
[l — on) 

where Fid is the free energy of the ideal gas: 

^ = -c„ In T + In n . 

The comparison of the calculation the binodal using (^1]) with 7 = (stan- 
dard incompressible volume) and 7 = 0.03 are presented in Fig. ^ Note 
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that the curve for Sd with 7 = was shifted for matching the characteris- 
tic minimum of these curves (see Fig. |^). The influence of the compressible 
proper molecular volume b on the binodal is shown in Fig. |^. As expected 
only the liquid branch (upper curve) is sensible to the dependence of b on 
the parameters of the state. The difference between the liquid branches of 
the binodal with 7 = and 7 = 0.03 increases with increase of the density. 
At t < 0.5 the role of the compressibility of the proper volume becomes very 
essential. For the model ( pTD the equation of state (EoS) as well as vapor- 
liquid phase equilibrium can be obtained. The results of the calculations are 
shown on Fig. |^ and display good correlation with the experimental data. 
The difference between the proposed model of compressible volume and the 
standard model with b = const is illustrated on Fig. |^ 

Similar results are obtained for other liquids with linear molecular geom- 
etry {i = 5) and low dipolar moments like N2, F2, NgO, COS for which the 
difference between low temperature values of the specific heat in coexistent 
phases does not exceed 2. 

The liquids with polyatomic molecules i = Q with hydrogen groups, like 
H2S, NH3 and CH4 also show satisfactory correspondence with the data. 
For their homologs where the hydrogen atoms are substituted by some other 
element, e.g. SO2, NF3, CF4 the model is less accurate in representing the 
behavior of Sd- It should be noted that in cases like NH3 and HgO one could 
expect the greater deviations because of the strong associative properties of 
these liquids The main source of error for such liquids comes from the 



absence of the proper density dependence of the heat capacity contribution 
due to dipole-dipole interactions which is clearly seen from the comparison 
of the behavior of the gas branch of the heat capacity (see Fig. pUj). 

Note that in cases of SO2, NF3, CF4 the parameter 7 takes the value 
7 ^ 0.025, which is smaller than that for their hydrated analogs H2S, NH3, 
CH4. It correctly represents the meaning of 7 as the compressibility of the 
free volume. Indeed, from the physical point of view it can be explained based 
on the fact that the hydrogen group has lesser moment of inertia and thus 
supports almost free rotation of the molecule in a cell, which result in greater 
compressibility. Therefore for the molecules like SO2, NF3, CF^ the influence 
of the neighbors is significant. For such liquids with polyatomic molecules 
further improvement can be reached with the help of the modification of the 
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particle-in-a-cell approach where the contribution of the molecular 

motion in a cell formed by the neighbors to the specific heat is taken into 
account (see Fig. |TTp. Also the molecular shape could be taken into account 
using the extension of the CS approximation for the nonspherical molecules 



4. Discussion 



In this work the main attention is focused on the behavior of the diameter 
Sd{t)of the vapor-liquid binodal in term of the temperature-entropy for dif- 
ferent molecular liquids. It is shown that Sd{t) for low molecular substances 
with non-spherical molecules has the surprising non-monotonous behavior. 
It is essentially different from that for the diameter ridit) of the binodal in the 
terms of the temperature- density. The latter is a monotonous function of 
the temperature excepting, maybe, the narrow vicinity of the critical point. 
Such a difference is connected with more complicated physical nature of the 
entropy. It is the sum of the two independent contributions: one of which 
is proportional to the density and the second is formed by the irreducible 
many particles correlations as well as the rotational and vibration degrees 
of freedom. At the same time the peculiarities of the diameter nd{t) of the 
binodal is mainly determined by the thermal expansion effects, which are 
very different for liquid and vapor phases. 

In the work it is shown that for low molecular liquids the non-monotonous 
behavior of Sd{t) is caused in the first place by the rotational motion of non- 
spherical molecules. The molecular motion influences also on such important 
molecular characteristics as the excluded volume in the equation of state for 
a system. Due to this the behavior of the entropy diameter is strongly con- 
nected with the EoS. The modified van der Waals EoS with the excluded 
volume, which is dependent on the external pressure is proposed. The cor- 
responding generalization of the ideal gas expression for the entropy is given 
too. It was shown that the diameter of the entropy and the vapor and liquid 
branches of the binodal are self-consistently and successfully described on 
the same basis. In connection with this we want to emphasize that the self- 
consistent reproduction of the PV^T-data and the diameter of the entropy is 
the weight criterion for the modeled EoS. 

Our consideration of these questions leads to necessity to construct the 
rigorous statistical theory for the description of the rotational motion of 
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molecules and, in particular, for the clear establishing of the P, T-dependence 
for the excluded volume in an EoS. 

The many-particles contributions to the entropy are especially important 
for liquid metals, first of all for Cs, Rb and Hg. Here the difference in 
the interparticle interaction is caused by the strong reconstruction of their 
electronic subsystem. Therefore the analysis of the temperature dependence 
of Sd{t) for them is also very important. 
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